Fermion Particle Production in Dynamical Casimir Effect in a Three Dimensional Box 
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In this paper we investigate the problem of fermion creation inside a three dimensional box. We 
present an appropriate wave function which satisfies the Dirac equation in this geometry with MIT 
bag model boundary condition. We consider walls of the box to have dynamic and introduce the time 
evolution of the quantized field by expanding it over the instantaneous basis. We explain how we can 
obtain the average number of particles created. In this regard we find the Bogliubove coefficients. 
We consider an oscillation and determine the coupling conditions between different modes that can 
be satisfied depending on the cavitys spectrum. Assuming the parametric resonance case we obtain 
an expression for the mean number of created fermions in each mode of an oscillation and their 
dynamical Casimir energy. 
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I. INTRODUCTION 

The Casimir effect is regarded as one of the most strik- 
ing manifestation of vacuum fluctuations in quantum 
field theory. The presence of reflecting boundaries al- 
ters the zero-point modes of a quantized field, and results 
in the shifts in the vacuum expectation values of quan- 
tities quadratic in the field, such as the energy density 
and stresses. In particular, vacuum forces arise acting 
on constraining boundaries. The particular features of 
these forces depend on the nature of the quantum field, 
the type of spacetime manifold and its dimensionality, 
the boundary geometries and the specific boundary con- 
ditions imposed on the field. Since the original work by 
Casimir in 1948 [l[ many theoretical and experimental 
works have been done on this problem (see, e.g., and 
references therein). A new phenomenon, a quantum cre- 
ation of particles (the dynamical Casimir effect) occurs 
when the geometry of the system varies in time. Theoret- 
ical investigations in this regard are for various fields, ge- 
ometries, boundary conditions, and different dimensions. 
Here we mention just some of them: motion of a sin- 
gle reflecting boundary [(|, the vacuum stress induced 
by uniform acceleration of a perfectly reflecting plane 
[7|, a sphere expanding in the four-dimensional space- 
time with constant acceleration investigated by Frolov 
and Serebriany [H, 0_in the perfectly reflecting case and 
by Frolov and Singh [101 ] for semi-transparent boundaries, 
and more general cases of motion by vibrating cavities 
considered on the base of various perturbation methods 
[IiT - [Tc| . Particle creation from the quantum scalar vac- 
uum by expanding or contracting spherical shell with 
Dirichlet boundary conditions is considered in [l9j . In 
another paper the case is considered when the sphere ra- 
dius performs oscillation with a small amplitude and the 
expression are derived for the number of created particles 
to the first order of the perturbation theory [2(| ■ Consid- 
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ering that Creation of particles and dynamical Casimir 
energies of configurations depend on the nature of the 
particular quantum field, now in the present paper by 
using the result of 21( we concentrate on the creation of 
particles for the Dirac field. Using the Dirac field in three 
dimensional box we investigate the number of fermion 
production in this geometry with dynamical boundaries 
for an arbitrary motion of the walls and then for an os- 
cillatory modulation. 

In order to obtain the number of produced fermions or- 
ganization of this paper is as follows, in sec. II. we present 
an appropriate wave function which satisfies the Dirac 
equation in side a box with MIT bag model boundary 
condition. We consider the walls of the box to have dy- 
namic, then we introduce the time evolution of the quan- 
tized field by expanding it over the instantaneous basis. 
Following the steps given in (l8j we arrive at an infinite 
set of coupled differential equations for the coefficients of 
the expansion. We explain how we can obtain the average 
number of particles created after the end of the motions. 
Dependence of this number to the Bogliubove coefficients 
leads to find the Bogliubove coefficients in sec. III. We 
consider an oscillation with small amplitudes of oscilla- 
tions and determine the coupling conditions between dif- 
ferent modes that can be satisfied depending on the cav- 
itys spectrum. Assuming the parametric resonance case 
we obtain an expression for the mean number of created 
fermions in each mode of an oscillation and dynamical 
casimir energy corresponding to the created particles. 



II. EXPANSION OF DIRAC FIELD OVER THE 
INSTANTANEOUS BASIS 

Let's have a brief review on the construction of the 
eigenfunctions of the Hamiltonian for the Dirac field in- 
side a box [2l| . In order to satisfy the purpose of appro- 
priate confinement physically relevant boundary condi- 
tion must be imposed. A proper boundary condition for 
the Dirac field is the MIT bag model boundary condition. 
This model was first considered by Bogolioubov [22[ and 
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later developed as the MIT bag model by Chodos ct al. 
[23| for hadrons. It is usually said to imply that there 
is no flux of fcrmions through the boundary. However, 
it implies an even stronger condition that is the absolute 
confinement of the Dirac field ( see also [24[ ). Consid- 
ering all of the back and forth terms inside the box the 
most general stationary solution is 



solution introduced in Eq. ([T]). When the modulation 
begins the boundary conditions on the moving walls be- 
come time-dependent. To satisfy these time-dependent 
boundary condition we expand the mode functions with 
respect to an instantaneous basis (25| : 



u n (x, t) = ^2(h u d v + m)Q[ n) (% fe (x, f) 



(6) 



$(xi,X2,x 3 ,t) = i)(x 1: X2,x 3 )exp(-iEt), 
where the spatial part is 

a i(pixi+p 2 x 2 +p 3 x 3 ) 
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We suppose that two-component spinors are identical for 
simplicity and consequently we have an overall normal- 
ization coefficient /, 
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where 
(2) 
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U+ = T], li_ = 



-Ti- 



en 
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Here p denotes momentum operator and a, f3, rj, (i, v, 
t, and p are general two-component spinors and co- 
efficients / through r can be determined. Imposing the 
prevalent form of the MIT bag model boundary condition 
on the dirac field inside a cubic box of side a as follow 



Considering continuity of each field mode and its time 
derivative at t = the initial conditions are 



(1 ± ZT^VKX) 



= 0, i =1,2,3 



(3) 



= ±a/2 



Q k n \o) = -iw n f 5 n , k , 
and due to the normalization condition 12611 we have 



(9) 



on the Eq. ^ yields in the quantization condition for all 
components of the momentum 



/ = 



(E + m) 
2EV '' 



Pi cotpia = —m, i = l,2, 3. 



(4) 



Now suppose that the distance a varies as a function of 
time a(t). Assume that a(t) has a constant initial value a 
and after a time AT the modulation stops and a(t) takes 
its initial value. The Fourier expansion of the field for 
an arbitrary moment of time, in terms of creation and 
annihilation operators, can be written as 



where V denotes volume of the box. The expansion in 
Eq. © for the field modes must be a solution of the wave 
equation. Inserting it in the Dirac equation and taking 
into account that the (f k 's form a complete and orthog- 
onal set of solutions of the wave equation and that they 
depend on t only through a(t), we obtain a set of coupled 
equations for Q k n \t): 



V;(x, f) = ]T aj>„(x, t) + 6^„(x, t). 



(5) 



We use the index " in" to mean times before modulation 
of the system and a™ and b% n are the annihilation and 
creation operators correspond to the particles in the " 
in" region. The mode functions u n (x,t) form a complete 
orthonormal set of solutions of the wave equation with 
MIT bag model boundary conditions. For the static cav- 
ity each field mode is determined by the Eq. (0| and 
the mode functions u n (x,t) have the form of the static 



with u>l 



III. 



Q[ n \t)M^t) +u k 2 (t)Q ( k n) (t)ip k (x,t) = (10) 
2 ^'fc| 2 . 



m 



DERIVATION OF BOGOLIUBOVE 
COEFFICIENTS 



When the boundaries of the box return to their initial 
position (p k s are time independent and the right-hand 
side in Eq. tfTU]) vanishes. We call this region of time " 
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out" region and wc can define a new set of creation and 
annihilation operators. The Fourier expansion of the held 
after modulation is 



^(x,i) 



Y,< Ut Mx,t)+bl out v n (x,t). (11) 



In this region Eq.(fT0|) reduces to the following simple 
form 



where the summation index s runs over the spin states. 
We are not able to sum over quantized mode of a mas- 
sive Dirac field given by Eq. ((4]) analytically. By setting 
m = we concentrate on the massless Dirac field, the 
quantization condition Eq. Q in this case yields the sim- 
ple form of pi = (in + 5)— for a mode p and then conse- 
quently the dynamical Casimir energy for these massless 
case becomes 



(d t 2 +oj p 2 )Q p n \t) = 0, 
and then the solution reads 



)(«) 



(t) = 4 n) e^' + B' n ) f 



(12) 



(13) 



with A p n ^ and B p n ^ constant coefficients to be deter- 
mined. The creation and annihilation operators for par- 
ticles and anti-particles in the in and out regions obey 
the usual anticommutation relations. Using the Bogoli- 
ubov canonical transformation one can expand the out 
operators in terms of the in operators 



a k af + f3 k b 



■ in 

-k 

-l%afc + cf k b% 



(14) 



Substituting Eq. ([13]) in Eq. © and then in Eq. (fTT|) . 
and considering Eq. (|14[) we expand ipk for out region in 
term of the creation and annihilation operators in the in 
region. After some calculation and by means of following 
relations that can be seen easily 



(<? • p)U- = (E — m)u+ 
(<? ■ p)u + = (E + m)u- 



eventually we get 



mu + + uj p u- - (E - m)u + \ (n) 
ma_ + u! p u + + (E + m)u- I p 



(15) 



(16) 



Using the relation between Bogoliubov coefficient for 
fermionic fields 



+00 



E = E E 



(20) 



s,n ni ,ri2 ,713— 



where u p = yj (m + i) 2 + (n 2 + i) 2 + (n 3 + \) 2 . 

Up to this point the equations are valid for an arbi- 
trary motion of the boundaries of the box. We only 
assume a(0) = a. We arc interested in the number of 
fcrmions created inside the box, so we look for harmonic 
oscillations of the walls which could enhance that number 
by means of resonance effects for some specific external 
frequencies. So we study the following oscillations 



a(t) = a(l + esin(fit)). 



(21) 



For small amplitudes of oscillations e <C 1 contenting 
ourselves with first order, the equations for the modes 
Ea.lflO]) takes the form (27 1 



)W(i) + W /QW(i) 



E + m 
2EV 



2eilcos(m)Y,9pkQr(t) 



-eft 2 sin(Ot) ^ g pk Q { k n) (t) + 2euj p 2 ma(SU)Q^ (t) 



(22) 



where 



1, 



(17) 



the mean number of particles produced in the mode p 
through an arbitrary modulation of the single fermion 
mode is the average value of the number operator with 
respect to the initial vacuum state 



N„ 



l&l 



l-2u 2 p \B^\ 2 . 



(18) 



In the presence of the boundaries, all of the components 
of the momentum are subject to quantization condition 
Eq.Q. Therefore the dynamical Casimir energy related 
to the particles production is given by 



E = EE n pV* 



\p\ 2 , 



(19) 



9 P k = a(t) / d"x 
Jq 



»(*) , s 3v>;(x,t) 



M*,t)- (23) 



Since s <C 1 it is natural to assume that the solution 
of Eq. is of the form 



Q(»)(t) = A p n \t) e^ 4 + B^\t) 



(24) 



with the function A p (t) and B p L> (t) varying slowly with 



time. We insert Eq. 



into Eq. (|2"2"j) to obtain differ- 



ential equations for A p n \t) and B p n \i). After neglect- 
ing their second derivatives and multiplying equation by 
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we average over the fast oscillations, we have 



initial condition mentioned in Eq. (|9|), it reads 



dA 



in) 



dr 



E + m 
2EV 



+E 



wfc + 0) 
(25) 



(n) 



r/r 



£7 + m 
2EV 



A- 



/,■ 



N 



(26) 



E 

k 



where r = et is a time scale. Now we shall solve equa- 
tions ([25]) and (|26| . Depending on the walls frequency 
and the spectrum of the static cavity we have different 
kinds of solutions. If we consider VL = 2lo p namely the 
frequency of the boundaries are twice the frequency of 
some unperturbed mode, in this condition if u>k — lu p = CI 
the resonant mode p will be coupled to some other mode 
k. Let us suppose that coupling conditions |o;p±ajfe| = fl 
is not fulfilled. In this case and for a massless field the 
equations ([231) arL( i (EHJ) reduces to 



dA™ 




dr 


w p 


dB p n) 




dr 


AV p 



(27) 
(28) 

The solutions of these coupled equations must satisfy the 



A (n) = _^ sinh( i^ r) 
P y/ZV AV h 

B {n) = i^ cosh( ^ r) _ 



(29) 
(30) 



By using Eq. (|18[) the average number of produced 
fcrmions in the mode p is 



iV p = l-^cosh 2 (^r), 



(31) 



and dynamical Casimir energy of the created fcrmions is 
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v w 

IV. CONCLUSION 



E = 2 E (i-^«** a (^))«P- 

ni,n 2 ,n3=0 



In this paper we have discussed the particle creation 
for a Dirac field in a three dimensional box with the MIT 
bag model boundary condition. We have considered all 
boundaries of the box to modulate during a finite time 
interval. We have used the Bogoluibove coefficients in 
order to obtain the number of fermion created during 
the motion. It is worth mentioning that we derived the 
Bogoluiobovc coefficients to be tetrad for the Dirac field. 
We have taken into account the usual parametric reso- 
nance case (CI = 2u) p ). We have also computed the dy- 
namical Casimir energy for this case for the mentioned 
modulation and the parametric resonance case. 
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